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ASYMPTOTICS OF DIRICHLET EIGENVALUES AND 
EIGENFUNCTIONS OF THE LAPLACIAN ON THIN 

DOMAINS IN R d 

DENIS BORISOV AND PEDRO FREITAS 



Abstract. Wc consider the Laplace operator with Dirichlct bound- 
ary conditions on a domain in R d and study the effect that per- 
forming a scaling in one direction has on the eigenvalues and cor- 
responding eigenfunctions as a function of the scaling parameter 
around zero. This generalizes our previous results in two dimen- 
sions and, as in that case, allows us to obtain an approximation for 



, Dirichlet eigenvalues for a large class of domains, under very mild 



assumptions. As an application, we derive a three-term asymptotic 
expansion for the first eigenvalue of d— dimensional ellipsoids. 

1. Introduction 

In his 1967 paper [J] Joseph studied families of domains indexed by 
psj ! one parameter to obtain perturbation formulae approximating eigen- 

values in a neighbourhood of a given domain. Within this context, 
he derived an elegant expression for the first eigenvalue of ellipses 
q ■ parametrized by their eccentricity e, namely, 

(L1) \ux 

■ where Ai = Ai(0) is the first eigenvalue of the disk - to obtain the 

eigenvalue of ellipses of, say, area ir, for instance, this should be divided 
by y/1 — e 2 and Ai(0) be the corresponding value for the disk. The 



2 (3-^) e 6 + £>(e 8 ), as e -> 0, 
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coefficient of order e 6 in Joseph's paper is actually incorrect - we are 
indebted to M. Ashbaugh for pointing this out to us, and also for 
mentioning Henry's book [H] where this has been corrected. Although 
in principle quite general, the approach used by Joseph yields formulae 
which, in the case of domain perturbations, will allow us to obtain 
explicit asymptotic expansions only in very special cases such as that of 
ellipses above. The failure to obtain these expressions may be the case 
even when the eigenvalues and eigenfunctions of the original domain are 
known, as this does not necessarily mean that the coefficients appearing 
in the expansion may be computed in closed form. An example of this 
is the perturbation of a rectangle into a parallelogram, which Joseph 
considered as an example of what he called "pure shear." 

With the purpose of obtaining approximations that can be computed 
explicitly, in a previous paper we considered instead the scaling of a 
given two dimensional domain in one direction and studied the result- 
ing singular perturbation as the domain approached a segment in the 
limit [BFj . This approach may, of course, have the disadvantage that 
we might now be starting too far from the original domain. However, it 
allows for the explicit derivation of the coefficients in the expansion in 
terms of the functions defining the boundary of the domain. As was to 
be expected, and can be seen from the examples given in that paper, 
these four-term approximations are quite accurate close to the thin 
limit. A more interesting feature of this approach is that in some cases 
it also allows us to approximate eigenvalues quite well away from this 
limit, as may be seen from the following examples. The application of 
our formula to the ellipses considered above yields 

(1.2) \ 1 {e) = ^ + - + l +{- + -)E + 0(e*), ass-+0, 

where we now considered ellipses of radii 1 and e, e being the stretch 
factor. The error in the approximation is comparable to that in Joseph's 
formula, except that equation (11.21) is more accurate closer to the thin 
limit while ( ll.lj) provides better approximations near the circle. This is 
also an advantage, since it is natural for numerical methods to perform 
better away from the thin limit, but to have more difficulties the closer 
they are to the singular case, suggesting that our formulae may also be 
useful for checking numerical methods close to the limit case. 
As another application we mention the case of the lemniscate 
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for which we have 

w . 2tt 2 2y/3n 97 / 593 V3n\ 2 

A l(£ ) = - + — + ^ + + — J e+0(e% as s - +0, 

yielding an error at e equal to one which is in fact smaller than in the 
case of the disk above. For details, see [BFj . 

In the present paper we extend the results in |BF] to general dimen- 
sion, in the sense that we now consider domains in M d which are being 
scaled in one direction and approach a (d — 1) -dimensional set in the 
limit as the stretch parameter goes to zero. Due to the increase in 
complexity in the corresponding formulae as a consequence of the fact 
that we are now considering arbitrary dimensions, we only obtained 
the first three non-zero coefficients in the asymptotic expansion of the 
principal eigenvalue. However, because of smoothness assumption near 
the point of global maximum, these include the coefficients of the two 
unbounded terms plus the constant term in the expansion - we know 
from the two-dimensional case that lack of smoothness at the point of 
maximum will yield other intermediate powers of e [Fj IFSj . 

As an example, we obtain an expansion for the first eigenvalue of the 
general d— dimensional ellipsoid 

\ 2 




+ (-] <1 

a d , 

where the a^s are positive real numbers. If we choose as projecting 
hyperplane that which is orthogonal to the Xd axis we obtain 

2 1 

d i=i * 



;i.3) 



d— 1 - - d— 1 d— 1 
^ l=\ 1 1=1 ]=l+l J . 

+0(e 1 / 2 ), as e -> +0. 



Besides the added complexity of the formulae, there are now extra 
technical difficulties related to the fact that there may exist multiple 
eigenvalues requiring a more careful approach. As in the two dimen- 
sional case, the asymptotic expansions obtained depend on what hap- 
pens locally at the point of global maximum width. Also as in that case, 
we cannot exclude the existence of a tail term approaching zero faster 
than any power of e. However, we conjecture that if the boundary of 
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the domain is analytic, then the expansions will actually correspond to 
the series developments of the corresponding eigenvalues. 

In the next section we establish the notation and state the main 
results of the paper, which are then proved in Sections [3] and HI In the 
last section, and as an application, we derive the above expression for 
the first eigenvalue of the ellipsoid. 

2. Statement of results 

Let x = (x',Xd), x' = (xi, . . . ,Xd-i) be Cartesian coordinates in M. d 
and respectively, d ^ 2, and uj C M d_1 be a bounded domain 

having C 1 -boundary. By h± = h±(x') G C(uJ) we denote two arbitrary 
functions such that H(x') := h + (x') + h_(x') ^ for x' G uj. We 
consider the thin domain defined by 

Q £ := {x : —eh-{x') < Xd < £h + (x'),x' G uj}, 

where e is a small positive parameter. We assume that the function 
H(x') attains its global maximum at a single point x E uj and that there 
exists a ball B' s (x) := {x' : \x' - x\ < 5} such that h± G C°°(B' s (x)). 
Let Hq := H(x) and the Taylor expansions for H and h_ at x read as 
follows 

oo oo 

(2.1) H{x') =H + Y,Hi(x' -x), h_(x') = h + ^T J h l (x' -x), 

i=2k i=l 

where Hi and hi are homogeneous polynomials of order i, H2k(x' — x) < 
for x 1 x, and k ^ 1. 

Our purpose is to study the asymptotic behaviour of the eigenvalues 
and eigenfunctions of the Dirichlet Laplacian — in Q e . Let x = 
x( x ') ^ C 00 (M d_1 ) be a non-negative cut-off function equalling one as 
\x' — x\ < 5/3 and vanishing for \x' — x\ > 8/2. Denote fl s e := fl e PI {x : 
\x' — x\ < 5}. 

Let 



<n ■ — t' t T 3 



be an operator in L 2 (M. d ~ 1 ). The spectrum of this operator consists of 
countably many isolated eigenvalues of finite multiplicity having only 
one accumulation point at infinity [GJ Ch. IV, Sec. 46, Th. 1]. By 

A n ,i < A nj 2 ^ A nj 3 . . . 

we denote the eigenvalues of this operator arranged in non-decreasing 
order and taking the multiplicities into account. Denote by ^ n , m the 
associated eigenfunctions orthonormalized in L 2 (M. d ~ 1 ). It follows from 
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[Gj Ch. V, Sec. 43, Th. 2] that the functions ^ n ,m decay exponentially 
at infinity. 

Our main results are the following. First, we obtain a two-parameter 
description for the eigenvalues. 

Theorem 1. Let A = A n>M = A n ^ M+ i = . . . = A n>M+N _i be a N- 
multiple eigenvalue of G n for a given n G N. Then there exist eigen- 
values \ n ,m{.£) of — Aq , m = M, . . . , M + N — 1 taken counting mul- 
tiplicities whose asymptotics as e — > +0 read as follows 

OO 

(2.2) X n , m (e)=e- 2 c^ m) +e-"J2^' m) V J , V-=e a , a := 

(9 o\ J n H _ 71-2712 J™>™) _ . 

^0 

and —c^T} are the eigenvalues of the matrix with the entries 

2>n 2 n 2 H^(H 2k+1 ^ n ^ ttn.OiaCR^i), m, / = M, . . . , M + N - 1. 

The remaining coefficients are determined by Lemmas 1 3. 6i and 3. 7 . 

As in |BF] , for sufficiently small e this allows us to derive the asymp- 
totics for specific eigenvalues, and we give the explicit expansion for the 
first eigenvalue in terms of the functions H and /i_ in the case where 
H 2 is negative for x' ^ x. 

Theorem 2. For any N ^ 1 there exists Eq = £o(N) such that for 
e ^ 6q the first N eigenvalues of — A® are Xi tm {e), m = 1, . . . , N. If 



(2.4) k = l, i7 2 (x') = -~5>^, 

i=i 

the lowest eigenvalue Xn(s) has the asymptotic expansion 



(2.5) A u (e) = + — + c i M) + O^e 1 ' 2 
e z e 



■2 



d-l 



j=1 H ' 



4 M) = JfA( 3H i(0 ~ 2 ^4(0)*o^o) i2(Rd - 
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d-1 q\/4 



[2.6) *o(0:=n^ e "^ 



/ ±i Q^fl .dr. 

_ J PP3^3 

\P,3= 



m') ■= Me) E 



; 2H*6 j (2B p + 8 j 



d-i 



- E 



71" Ppqj^p^q^j 



J 3 ' 

p,q,j= 



■ sX Hl% + q + j ) i 

where it is assumed that H 3 (x') is written as 

d-i 

Hs(x ) ^ ^ fipqj£,p£,q^j j 

and £/ie constants f3 pq j are invariant under each permutation of the 
indices p, q, j : 

(2.7) (3pqj = (5pj q = Pgpj = (3 q j p = (3j pq = (3jqp- 

Remark 2.1. The assumption (12.41) for H 2 is not a restriction, since 
we can always achieve such form for H 2 by an appropriate change of 
variables. 

3. Proof of Theorem CD 

In this section we construct the asymptotics for the eigenvalues and 
the eigenf unctions of the operator — A^ £ . This is first done formally, 
and justified rigorously afterwards. In the formal construction we em- 
ploy the same approach as was used in \BF\ Sec. 3]. 

We are going to construct formally the asymptotic expansions for 
the eigenvalues A„ im (e), m = M, . . . , M + N — 1 which we relabel as 
\i m \ m = 1, . . . , N. We denote the associated eigenfunctions by tpe 
We construct their asymptotic expansions as the series 



oo 

i 



A M = £ - 2/i ( m) j = C M + J2 ct\ 



(3.1) 



i=2k 

oo 



i=0 



(d c\ d x' — x x d + eh_(x') 

77 sH(x') 

We postulate the functions ^| m ^(£) to be exponentially decaying as 
£' — > +oo. It means that they are exponentially small outside fif 
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(with respect to e) . In terms of the variables £ the domain fl s e becomes 
{£ : < ^~\0 < £ d < 1}. As rj -»• 0, it "tends" to the layer 
II := {£ : < £d < 1} and this is why we shall construct the functions 
ipi as defined on II. 

We rewrite the eigenvalue equation for ip e and A e in the variables £, 
(3.2) 



d-l 

£ 

i=i 



2fe+l 



d d d d 



d-l 



+ T) 



2k+2 



E 9 r^2 9 , „2k+2 u 



^ m) =0 on dU 



= /4 m ty< m) in n, 



where K x = Ki(£, rj), % = 0, . . . ,d, 

1 



H 2 (x + ri£ r 
1 

1 



dh- , ... „ OH . ... 

OX; OX; 



K Q (Z,ri) = -H-\x + r£ , )A x ,H(x + rit') 



Remark 3.1. We have introduced the factor y/H(x') in the series ( 13.11) 

for ipe in order to have a symmetric differential operator in the equa- 
tion (P - 

We expand the functions i^j into the Taylor series w.r.t. 77 and 
employ (12.11) to obtain 



(3.3) 



j=2k 



j=0 



Kf(0 :=Pf (0 + ^(0. 

00 

^o(£,r/) = E^ (0) (0, 



i=0 



d-l. 
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where Pf \ Q$ are polynomials, and, in particular, 



(3.4) 



r>(A(ti\ - 2 -^2fc(^ ; ) p (d) _ 2H 2 k+i{.j') 



We substitute (13.11) . (13.31) . (13.41) into (13.21) and equate the coefficients 
of like powers of t]. This leads us to the following boundary value 
problems for ip^ 



1 d 2 



(3.5) \Htde~ 



( m ) i „/,( m ) 



V>v" =0 in n, 

d 

ip( m) = on dU, j = 0, . . . , 2k - 1, 
1 <9 2 



(3.6) ( a gggfc(eO <9 2 (m) \ ( m ) . 



4r } = o on an, 



(3.7) 



/ 1 d 2 ( m )\ ,(m) f a 2H 2 k(£') d 2 ( m) \ (m) 

Use + c ° ) *> = \ A " - + C2t l 

j-2k-l 

+ cfVi m) + £ c£^J m) + *f° in n, 



£ 

^ (m) = on an, j ^ 2A; + 1, 

j'-2fc-l 
9=0 



d-1 3-2 



f 2^ K i— 2^ + ^'+2^2 + ^-2, 



,(0) 

at/ Xa iyj - s ~ 2 a^ ^ 1 j +2k f)f 2 ^ 1 

i=l s =0 

where = 0. Problems (I3.5P can be solved explicitly with 
(3.9) ^ m) (0 = *; m) (nsinvm6, c 



rr2 ' 
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where j — 0, . . . , 2k — 1, and ^y 1 ' are the functions to be determined. 

The last identity proves formula (12. 3p for c n . 

We consider the problem (13. 6p as posed on the interval (0, 1) and 
depending on It is solvable, if and only if the right-hand side is 
orthogonal to sin7m^ in £2(0, 1). It implies the equation 

(3.10) - (A, + frM = 4^)^M i n R^-i. 

Thus, is an eigenvalue of the operator G n , i.e., = A. Then 
is one of the eigenfunctions associated with A. These eigenfunctions 
are assumed to be orthonormalized in L 2 (M. d ~ 1 ). We substitute the 
equation (13 . 1 0[) into (13.61) and see that the formula (I3.9P is valid also 
for j = 2k. 

The problems ( 13. 71) are solvable, if and only if the right-hand sides 
are orthogonal to sinirn^ in L 2 (0, 1). It gives rise to the equations 

j-2k-l 

(3.11) (G n -A)¥^ 2k = f^+ c^y^ + c { ^^\ 

1 



*H = *J m, (0 := 2 1 ^ m) (0sinvr<,d^, 



(3.12) 







To solve the problems (13.71) . (13.111) we need some auxiliary lemmata. 
The first of these follows from standard results in the spectral theory 
of self-adjoint operators. 

Lemma 3.2. Let f E L 2 (R d ~ 1 ). The equation 
(3.13) (G n -A)u = f 

is solvable, if and only if 



(/,^ m) )L 2 ^-i) = 0, m = l,...,N. 

(m) 



The solution is unique up to a linear combination of the functions \l/ 
By g n we denote the sesquilinear form associated with G n , 

Q n [u,v] = (Vu, Vt)) ij(H d-i) - (H 2 kU,v) L2{R d-i). 

The domain of this form is 

V( Qn ) = W 2 1 (R d ~ 1 ) n {it \ (1 + \i'\ k )u G L 2 (M d ~ 1 )}. 
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By T>(G n ) we denote the domain of G n . The set C^°(M. d 1 ) is dense in 
V(g n ) in the topology induced by g n ([D] Ths. 1.8.1, 1.8.2]). 

Lemma 3.3. Let f G L 2 (M d_1 ), u G L 2 (K d_1 ) n W 2 1 (5 r ) /or eac/j 
bounded domain S C and for each <fi G T>(q) the identity 

(3.i4) y v M -v0d£'- y (H a (o-A)^dc= y ^df 

E d-i K d-i R d-i 

/ioWs irue. TTien it G T>(G n ) and the equation A3.13\) is valid. 

Proof. Let Xi = Xi{t) be a non-negative infinitely differentiable cut-off 
function taking values in [0, 1], equalling one as t < 1, and vanishing as 
t > 2. It is clear that for each t > the function w(f )xi(|f |£) belongs 
to £>(fl„). We substitute 0(f) = u(f )xi(|f |£) into 6&D and integrate 
by parts, 

IIXI^HLOR'*- 1 ) ~ ( H 2kXlU, Xl«)i 2 (Rrf-i) = (Xl/,Xl^)L 2 (IRd-i) 

(3.15) 
Hence, 



+ j ( uA £'Xi> M) L2(Kd -i) + AUxiullia^d-i). 

IIVtiHiafB/ (o)) - {H 2 kU,u)L 2 (B'(0)) <: 1 1 / 1 1 ^(R^- 1 ) 1 1 « I U 2 1 ) 

( 3 - 16 ) * „ +2|l „ ... I|2 

+ IPlUa^- 1 ) + A ll -u llL 2 (M d - 1 )' 

where the constant C is independent of t, and B' r (a) :={£': \£' — a\ < 
r}. Passing to the limit as t — > +0, we conclude that u G T>(g n ) and in 
view of (I3.14p this function belongs to T>(g n ) and solves the equation 
(JSH. □ 

Let 23 be the space of the functions / G C 00 ^!^" 1 ) such that 

8 T f 

for each r G Z+, 7 G Z H 



7d 



21 



Lemma 3.4. Lei / G 23, and u be a solution to A3.13\) . Then u G 23. 

Proo/. Since u G £>(G„), we have Vw G L^" 1 ), (1 + |6'| fc )w e 
L 2 (M d ~ 1 ), and due to standard smoothness improving theorems u G 
C°°(R<f-i). The identity 1ETI5]) is also valid with X i replaced by Xi{\£'\t)\£' f. 
Employing this identity and proceeding as in f)3.16p . we check that 
(1 + |£T)Vu g L 2 (M d_1 ), (1 + |f \ k+p )u G L 2 (R d - 1 ), if (1 + \?f)u G 
L 2 (M d ~ 1 ) for some /? G Z+. Applying this fact by induction and using 
that (l + |f \ k )u G L 2 (M d_1 ), we conclude that (l + |f HVw G L 2 (M d_1 ), 
(1 + |f \ k+ ^)u G L 2 (R d ~ 1 ) for each 7 G Z+. 
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We differentiate the equation (13.1 3D w.r.t. 

(G - -^wrw, + -air"- 

The right hand side belongs to L 2 (IR d_1 ) and the function J|- satisfies 
the hypothesis of Lemma 13.31 Applying this lemma, we see that G 
T>(G n ). Proceeding as above, one can make sure that 

Oil 

(l + imv^e^fR") 

for each 7 6 Z+. Repeating the described process, we complete the 
proof. □ 



As it follows from Lemma \3. 21 the solvability condition of the equa- 
tion fl3~TT|) is 

j-2k-l 
9=1 

where m, / = 1, . . . , N, and 5 m i is the Kronecker delta. Here we have 
supposed that the functions ^q"^ are orthonormalized in L 2 (lR d_1 ). In 
view of (I3.12p these identities can be rewritten as 

j— 2fc— 1 

(3.17) 2(if" ) ,4 /) ) i2(n) + 2 £ c^f ) ,rf ) )x 2 (n) + 4 m) ^ = 0, 

9 =1 

where m, I = 1, . . . , N. 

Consider the problem (13 .7p for j = 2k+l. The solvability condition 

is the equation (13. lip for the same j. Since \&q"^ G 53, the same is true 

f° r /Ife+i- By (13.171) . this equation is solvable, if and only if 

(3.18) T,2 fc+1) + c^Srm = 0, m, I = 1, . . . , N, 

The definition of C\ and (13.41) yield 

(3.19) Ti? +1) = 27rVtf - 3 (iW^ m) , ^X^y 

Hence, the matrix y( 2fc+1 ) with the entries T^ +1 ^ is symmetric. This 

matrix describes a quadratic form on the space spanned over \?o , 
m = 1, . . ., N. By the theorem on the simultaneous diagonalization 
of two quadratic forms we conclude that the eigenfunctions can 
be chosen as orthonormalized in L 2 (M> d ~ 1 ) and, in addition, so that 



L 2 (U) 
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the matrix x( 2fc+1 ) is diagonal. In what follows we assume that these 
functions are chosen in such a way. Then identities (13.181) imply 

(3-20) 4Si = -rr +1 \ 

where Tm ^ are the eigenvalues of x( 2fc+1 ). 

By Lemma 13.21 the solution to (13. lip for j = 2k + 1 reads as follows 

N 

(3.21) ^\0 = ®f n \0+J2tM\ 

P =i 

where <^ m) is orthogonal to all ^\ I = 1, . . . , N, in L 2 (M d ~ 1 ) and b$ 

are constants to be found. It follows from Lemma [3.31 that $i G 2J. 
The definition (13. 8p of C\ and the equation (13. lip for j = 2k + 1 imply 
that the right-hand side of the equation in (13.71) for j = 2k + 1 is zero. 
Hence, the solution to the problem (13.71) for j = 2k + 1 is given by the 

formula (EOty . where tfj^ is to be found. We substitute (j3^ . fl3T2TD 
into the equation (ETTTj) for j = 2& + 2. In view of fl3~T7l) and fl3~20l) 
the solvability condition for this equation is as follows 

,(m), (2*+l) _ (2*+1)n i M v 
22) C 2fc+2°m« 

+ 2(£ 2 4 m) + £ 1 $[ m) sin vr< d , rf)^^ = 0, Z = 1, . . . , iV. 

Assume that all the eigenvalues Tm k+l ^ are different. In this case the 
last identities imply 

. M _ 2(£ 2 4 m) +^^i m) sin 7 r< d ,4°) L2( n) . . 

(3.23) ri 2fc+1) - r/ 2fc+1) ' ™ # ' 

4^ = -2(£ 2 ^ m) + £M n) sinvr< d ,4 m) ) L2( n)> 

and we can also let b^™\ = 0. 

Now suppose that all the eigenvalues Tm ^ are equal. In this case 
the equations (13.221) do not allow us to determine the constants b\™ 
for m ^ I. Consider the matrix x( 2fc + 2 ) with the entries 

T^ +2) ■= 2(C 2 4 m) + ^ sinvma4°) i2(nr 
Lemma 3.5. The matrix x( 2fc + 2 ) is symmetric. 
Proof. Integrating by parts, we obtain 
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Since by (O) 

rAl) _ _ (}_^_ , >»A _ M ,(m) 

^1^0 — ^i/ 2 <9£ 2 J <t J 2k+l °2fe+ir0 ) 

in view of ([SI]), pi) . flQH) we have 

2($S m) sin 7 r< d ,£ 1 4 /) ) L2(n) = (Q™, (G n - A)*®)^^, 

= ((<?„ - A)$^, $«) L2(Rd _ 1} = 2(A4 m \ *<™> sin7r< d ) L2(n) . 

□ 

Since we supposed that all the eigenvalues of T^ 2k+1 ^ are equal, we 
can make orthogonal transformation in the space spanned over \&g m \ 
m = 1,...,N, without destroying the orthonormality in £2(11) and 
diagonalization of T^ 2k+1 \ We employ this freedom to diagonalize the 
matrix x( 2fc + 2 ) which is possible due to Lemma 13.51 After such di- 
agonalization we see that the coefficients c^ +2 are determined by the 
eigenvalues of the matrix T^ 2k+2 ^: 

S m ) _ __(2fe+2) 
u 2fe+2 'm 

If all these eigenvalues are distinct, we can determine the numbers b\ ^ 
at the next step by formulae similar to (13.231) . If all these eigenval- 
ues are identical, at the next step we should consider the next matrix 
rp{2k+3) anc j diagonalize it. 

There exists one more possibility. Namely, the matrix x( 2fc+1 ) can 
have different multiple eigenvalues. We do not treat this case here. 
The reason is that the formal construction of the asymptotics is rather 
complicated from the technical point of view and at the same time it 
does not require any new ideas in comparison with the cases discussed 
above. Thus, from now on, we consider two cases only. More precisely, 
in the first case we assume that the matrix T^ 2k+1 ^ has iV different 
eigenvalues Tm k+1 \ m — 1, . . . , N. In the second case we suppose that 
the matrix T^ 2k+1 ^ has only one eigenvalue r^ 2fc+1 ^ with multiplicity 
N, while the matrix y( 2fc + 2 ) has iV different eigenvalues Tm k+2 \ m = 
1,...,N. 

Lemma 3.6. Assume that the matrix T^ 2k+1 ^ has N different eigen- 
values and choose \Pq being orthonormalized in L 2 (M d ~ 1 ) and so that 
the matrix T^ 2k+1 ^ is diagonal. Then the problems Ii3.5\) , \3. 6\) , ^3. 7| ) 
have solutions 

N 

p=i 
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Here the functions ifjj are zero for j ^ 2k + 1, while for other j they 
solve the problems 



j-2k-l 

+ E c^f ) + ^ M -2(^ (m \sm7r< rf )L 2 (o,i)sinvr< d m n, 

<j=2fc+2 

^ (m) =0 on 011, 
and are represented as finite sums 



^nere ^ G 23, V&? 2 G C7g°[0, 1], V>S(0) = ^(1) = 0, and the 

functions ipj™^ are orthogonal to sin7rn£d in 1/2(0, 1). Tne functions 

are solutions to the equations A3.11\) and are orthogonal to 

all ^\ I = 1,...,N, in L 2 (M d_1 ). The constants and cj m) are 
determined by the formulae 

2fF (m) ?//' )n I 4- V r (m) h {m) 

J.I ~~ (m) _ (0 ' TIT, y£ I, J ^ i-, 

T 2k+l T 2k+1 
S m ) _ A r M _ _(2fc+l) 

c M = _ 2(j p ? M^M )L2(n)) ^2fc + 2, 

j-2fc-l j-2fc-2 AT 

g=0 g=0 p=l 

Lemma 3.7. Assume that all the eigenvalues of the matrix T^ 2k+1 ' are 
identical and that the matrix T^ 2k+2 ^ has N different eigenvalues, and 
choose \E f Q m ' ) being orthonormalized in L/2(M. d ~ 1 ) so that the matrices 
T( 2fc+1 ) and T ( - 2k+2 '> are diagonal. Then problems RTfy) . [3~R) . (3~1\) 
have solutions 

^(0 =^ (m) (0 + ^ m) (nsmvr< rf 

N N 

+ E b ^iA\o sin ™^ + E b ^o\o- 
p=i p=i 
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Here the functions ifjj 7 ^ are zero for j 2k + 1, while for other j they 
solve the problems 

1 d 2 ( m )\ ~( m ) / 2H 2 k(C') 9 2 \ ~( m ) 

j-2fc-l 

+ E cH^+^-S^H^inTrn^^^sinTrned in n, 

q=2fc+2 

^ (m) = o on an, 

3~2fc-l 

F M. = ^ ^.^(^M + vpMsinvrn^) 

Ar j-2k-2 

+ E E ^l P A^M p w< rf 

p=l q=l 
Ar j-2fc-3 

+ E E 

p=l q=0 

and are represented as finite sums 

^ m) (o = E^S(w£j,(&), 

^ere ^ e 2J, ^3 <= C °°[0,1], ^ 2 (0) = ^ 2 (1) = 0, and the 
functions ipj™ 2 are orthogonal to simm^d in L 2 (0, 1). The functions 
are solutions to the equations 

j-l j-3 N 

(G n - a)*™ =j%i + e +EE 4:U«*o p) 

g=l g=l p=l 

3-2 TV Af 

g=l p=l p=l 

and are orthogonal to all \&g , I — 1 : . . . : N, in L 2 (M d ~ 1 ). The constants 
bjp and Cj are determined by the formulae 



Z \ r j+2k+2i YQ ) + 2^ c j+2k-q+2 U q,l 

(m) _ (i) , m T h J ^ 1 > 

T 2fe+1 T 2fc+1 
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-M _ A r (m) - _^r( 2fc + 1 ) r (m) - _^( 2fc + 2 ) 



C 2fc ~~ A ' C 2k+l ~ T i C 2k+2 

C M = _ 2 (^ m) ,4 m) ) i2( n), j^2fc + 3. 

These lemmata can be proven by induction. 

Remark 3.8. We observe that if A is simple, then N = 1 and the 
hypothesis of Lemma 13.61 is obviously true. 

We denote 

^ j /M ( x ' ~ w x d + eh-{x') 



j=0 



7] EH(x' 



A M := 5" 2 c ( m) + ^ ^cf \ s > 2k. 

j=2k 

The next lemma follows from the construction of the functions ipj^ 
and the constants cj . 

Lemma 3.9. The functions tpi™ solve the boundary value problems 
(3.24) -(Ag,+AH)^ = #, m = l,...,N, 

where the right-hand sides satisfy the estimate 
(3-25) ||2i?|k(Q £ ) = 0(rj*- 3J ^-% m = l,...,N. 

We rewrite problem f!3.24p as 

1 + A e 



where A e := (— + l) -1 . This operator is self-adjoint, compact and 
satisfies the estimate \\A e \\ ^ 1. In view of this estimate and (13.251) we 
have 



-A o (m) 



1 + Ag? ^ 



<C m , sV s -^ +2k , m = l,...,iV, 



L 2 (n E ) 

where C m , s are constants. We apply Lemma 1.1 to conclude that there 
exists an eigenvalue {e) of A e such that 

\g<r\e) - (1 + AM)" 1 ! ^ C m ,^"^ +2fc , m = 1, . . . , N. 

Hence, the number \i rn \e) := (^"^(e)) 1 — 1 is an eigenvalue of the 
operator — A^ , which satisfies the identity 

(3.26) | A M( £ )-AM| ^C m ^ s - 7J ^-\ m = l,...,N, 
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where C m)S are constants. 

Let be a monotone sequence such that C m ^ s rj ^ C m;S _i as e ^ 



e s ■ We choose the eigenvalue A e := Xj^ as e G [£« , £g+i). The 



inequality (I3.26P implies that the eigenvalue Ae has the asymptotics 
(12.21) . We employ Lemma 1.1 in |OSYl Ch. HI, Sec. 1.1] once again 
with a = C m ^r] s ~, d = ^/a. It yields that there exists a linear 
combination ipi m \x,e) of the eigenfunctions of — associated with 
the eigenvalues lying in [Xi™^ — d, Ae + d] such that 



Since the functions ^4,s are linearly independent for different m, the 
same is true for ipi m \-,e), if s is large enough. Thus, the total multi- 
plicity of the eigenvalues Ae is at least N. The proof is complete. 



In order to prove Theorem [2] we need to ensure that, for sufficiently 
small e, the asymptotic expansions for Ai >m , m — 1, . . . , N provided by 
Theorem [1] do correspond to the first N eigenvalues of — (counting 
multiplicities). In [BF] this was done by means of adapting the proof 
of Theorem 1.1 in |FSj from the situation where h- — to our case. 
In the present context we need to show that, under the conditions 
for h±, this result may be extended to d dimensions. There are two 
important points that should be stressed here. On the one hand, we 
are assuming C°° regularity in a neighbourhood of the point of global 
maximum, and thus do not have to deal with what could now be more 
complex regularity issues at this point. On the other hand, since the 
proof of eigenvalue convergence given in |FSj is based on convergence 
in the norm, it is not affected by details related to the possible higher 
multiplicities as was the case in the derivation of the formulae in the 
previous section. 

While still using the notation defined in Section [2l we also refer to 
the notation in |FSj . In particular, the function h and the operator H 
defined defined there correspond to our width function H and operator 
G n , respectively. We begin by assuming H to be strictly positive in UJ. 
Let thus 



UTK-,e)-4^\\ LM = o( v 




4. Proof of Theorem [2] 



ip(x',x d ) = ^ x (x',x d ) = x(x') 



eH{x') 



2 



sm 



ix(x d + eh-{x')) 
eH(x') 
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As in IFSI. we have 



U X ( X ^ X d)\\L 2 (Q £ ) = I X 2 (x')dx' 



while now 

reh+lx') 



\ \ \V^ x (x' jXd )\ 2 = / |V X (:r')| 2 + 

J uj J — eh— (x') J uj 



TT 



- + v(x') x 2 {x')dx' 



e 2 H 2 (x 



with 



v[x 



TT 



H 2 (x') 



-VH(x') - Vh_(x') 



2 1 






+ 4 


3- 





VH(x 



/\|2 



In the notation of [FS] . the potential W e appearing in the quadratic 
form q e [x] (equation (1.4) on page 3 in that paper) is now defined by 



W £ {x') 



71 



2 r 



1 



H 2 (x') 



H 2 (x) 



+ v(x'). 



We consider the scaling x' = e a t as before, which causes the domain uj 
to be scaled to uj £ = e a u. Then the proofs of Lemma 2.1 and Theorem 
1.2 go through with minor changes (note that m = 2k, while / and I e 
should be changed by uj and u e , respectively). Similar remarks apply 
to the proofs in Section 4 of |FSj leading to the proof of Theorem 1.3, 
except that due to regularity we do not need to worry about separating 
the domain into different parts as was necessary there for the intervals 
h. 

Finally, we relax the condition on the strict positivity of H mentioned 
above. This again follows in a similar fashion to what was done in 
Section 6.1 of [FS]. 

We are now in conditions to proceed to the proof of (12.51) . In the 
case considered the lowest eigenvalue of G\ is A = Ylj=i ®h wn il e the 
associated eigenfunction is given by (12.61) . This proves the formula 
for «4 . In view of Remark 13.81 we can employ Lemma 13.61 to cal- 



culate cjj 1 , c^' 1 ^ 
and # 3 (-0 = - 
Theorem [1] it yields that c 3 



Since \l/o is even w.r.t. each £j, \ 
■H 3 (£'), we conclude by fl3~T9|) that T- 



1. 

(3) 
11 



,d - 1, 
0. By 



0. 



(4.1) 



The equation (13. lip for with j = 3 reads as follows 

2tt 2 



(Gi-A)tf a 



H 3 



H 3 V . 
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We seek the solution as \&i = P^O; where R is a polynomial of the 
form 

d-l d-l 

(4-2) R(g) := E + E <% 

where C P9J -, Cj are constants to be found, and C pg j are invariant under 
each permutation of the indices p, q, j. We also note that such a choice 
of R ensures that ^cOl^r 1 *- 1 ) = 0- We substitute (14. 2 p and the 
formula for ^ into (14. ip taking into account (12. 7p . 

d-l d-l d-l 

2 ^2 ( P + d q + d j) C PQj£p£q£j + 2 E ^/'M + 6 E 

p,g,j=i j=i p,i=i 

2tt 2 d ~ l 

^1 ftpqj£,p£,q£,j ■ 



It yields the formulae 



P.9J = 1 



1 j ra mVp + Oi + OiY ° J 2 1^0,(20, + e^' 

It is easy to check that 

,, 4) ^<0 = -£ff<f), ^ = ^ 

Pf (£') = V(3A|(0 - 2H H,(e)). 
Employing these identities, we write the formula for c^ 1 ' 1 ^ from Lemma l3~6l 
c (i,i) _ _ 2(F 4 ,^ ) L2(n) = -2(£2^o,^o)L2(n) - 2(£ 1 V 1 sin vr£ d , Vo)L 2 (n) 

r 2/DM,T,_ ,T, A . . .. /-d(0). 



=7T 2 (P 4 W ^0,^0)L 2 (R d -i) - (P W ^ ,^o 



- 1 ) 



+ 4vr E ( Q? sin ^o£d cos 7r£ d J 

i=i V a * / L 2 (n) 

- 2E « ») 2 (SM + - 2 (^f5i, 

i=i V / L 2 (n) 

- 2 (pf*o, + ^ E «. ? - E *< 

2H ° i=l i=l V d & 



20 



7T*((Pf +P 3 W )*0,*0 
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2 d-1 



/ 9hi 







i=l 



We substitute (JO)) , (lOj) into this identity and arrive at the 

desired formula for c 4 . 



5. The ^-dimensional ellipsoid 



As an application of our results, we will derive the expression (11. 3p for 
the asymptotic expansion for the first eigenvalue for a general ellipsoid. 
From the equation defining the boundary of £ and assuming that, as 
mentioned in the Introduction, we are doing the scaling along the Xd 
axis, we have 



h+(x') 



ad 



1-1^ 

a 1 



ad-i 



1/2 



while H(x') = 2h±(x'). We thus have x located at the origin and 
H = lad- Expanding H around x we have 



H(x') 



2a d - a d 



ad 
'T 



a\ ) T l a d -x 



El) + 
ax) ^ 



Xd-l 

a d -\ 



9 ( xxx 2 \ 2 , 9 / X1X3 V 



ad-2ad-i 



+ ..., 

yielding H k = hk = for odd A; and 



Hence 



and 



Ho(x') 



H,{x') 



d-1 , x 2 
Xa 



-a d J2 

8=1 



1 d-1 



t=l j=l 



0C %tC j 



a _ 7T 
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Note that since if 3 is identically zero, there is no need to compute 
It is now straightforward to obtain 

7T 2 7T 1 

Co = — o and c*' 1 ' = — . 

4a 2 2 2a d *-'a i 

It remains to compute 




which, after some further simplifications, yields the desired result. 
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